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Abstract 

This paper investigates the Arnowitt-Deser-Misner (hereafter ADM) form of spherically sym- 
metric gravity with variable Newton parameter G and cosmological term A c . The Newton param- 
eter is here treated as a dynamical variable, rather than being merely an external parameter as 
in previous work on closely related topics. The resulting Hamilton equations are obtained; inter- 
estingly, a static solution exists, that reduces to Schwarzschild geometry in the limit of constant 
G, describing a Newton parameter ruled by a nonlinear differential equation in the radial variable 
r. A remarkable limiting case is the one for which the Newton parameter obeys an almost linear 
growth law at large r. An exact solution for G as a function of r is also obtained in the case of 
vanishing cosmological constant. Some observational implications of these solutions are obtained 
and briefly discussed. 
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I. INTRODUCTION 



The last ten years have witnessed an encouraging progress in the application of 
renormalization-group methods to the nonperturbative renormalization of Quantum Ein- 
stein Gravity [l[ M [J 4]. In general, field theories which are nonperturbatively renormaliz- 
able are constructed by performing the limit of infinite ultraviolet cutoff at a nonGaussian 
renormalization group fixed point in the space of all dimensionless couplings g± which 
parametrize a general action functional. In the case of general relativity, an effective average 
action has been built [l| , and such a nonGaussian ultraviolet fixed point has been found in 
the case of the Einstein-Hilbert and higher-derivative truncations [2J. 

Several cosmological applications of this framework have also been considered. For ex- 



ample, in jl]] (see also 



it has been argued that the resulting scale dependence of the 



Newton parameter at large distances might mimic the presence of dark matter at galactic 



and cosmological scales. On the other hand, in early work by some of us [7] , we had tried to 
build an action functional where the running of the Newton parameter is ruled by suitable 
Euler-Lagrange or Hamilton equations, while being compatible with the renormalization- 
group flow in the neighbourhood of an ultraviolet fixed point. For this purpose, one adds 
to an action of the Einstein-Hilbert type (but with variable G, so that it is brought within 
the integrand) two compensating terms such that the action reduces to the York-Gibbons- 
Hawking 9] action for fixed G and A c , and takes the same functional form as in the ADM 
formalism for general relativity (see Eq. (2.1) below). 

The work in Ref. [7] focused on cosmological models with Friedmann-Lemaitre- 
Robert son- Walker symmetry, but of course other symmetries are also relevant in the in- 
vestigation of the early universe. In particular, we are here concerned with the requirement 
of spherical symmetry, for which a thorough Hamiltonian analysis in general relativity was 



performed, for example, in Ref. [10,]. In that case, the starting point is of course the 
Schwarzschild line element written in the curvature coordinates (T,R), i.e. 

ds 2 = —F(R)dT 2 + F-\R)dR 2 + R 2 {dd 2 + sin 2 6d<j) 2 ), (1.1) 

where, in c = 1 units, F(R) = 1 — 2 ^-- A space-time foliation is then introduced, with R 
and T smooth functions of new independent variables r and t. 

At this stage, section 2 generalizes this construction to models with variable Newton 
parameter G and cosmological term A c , obtaining the general Hamilton equations for such 



models and proving that a static solution, compatible with the fixed-point hypothesis, ac- 
tually exists. Section 3 studies the resulting nonlinear differential equation for G, while 
some phenomenological implications are analyzed in section 4, and concluding remarks are 
presented in section 5. 



II. ADM ACTION AND HAMILTON EQUATIONS WITH SPHERICAL SYM- 



METRY AND VARIABLE G, A c 

From the analysis in Ref. we know that the ADM action for a theory of pure gravity 
where the Newton parameter G and the cosmological term A c evolve in space-time as a 
result of renormalization-group equations in the early Universe, can be taken to have the 
form 



S 



16tt 



M 



K ij K ij -K 2 + ( VR-2A C 



-1*9 



NVhd 3 x dt, 



(2.1) 



G ™ G 3 

where Kij is the extrinsic curvature tensor of the spacelike hypersurfaces £ which foliate the 
space-time manifold, while is their scalar curvature, and \i is an arbitrary dimensionless 
parameter. In other words, it is possible to generalize the standard ADM Lagrangian and 
regard G as a dynamical field obeying and Euler-Lagrange equation, the underlying idea 
being that all fields occurring in the Lagrangian L should be ruled by L in the first place. 
This makes it possible to fully exploit the potentialities of the action principle. For this 
purpose, one adds to an action of the Einstein-Hilbert type (but with variable G, so that G 
is brought within the integrand) two compensating terms 7] such that the action reduces to 
the York-Gibbons-Hawking [s),!^] form for fixed G and A, and takes the same functional form 
as the ADM action for general relativity. Non- vanishing values of /i in (2.1) (cf [111 ]) ensure 
that no primary constraint of vanishing conjugate momentum to G arises. Its preservation 
in time (which is necessary because G is here a dynamical variable, on the same ground of 
the space-time metric) would lead, following the Dirac method [12], to further constraints 
with considerable technical complications. Such an action is sufficient for the purposes of 
an ADM analysis, which is what we do hereafter, but its generalization to a form invariant 
under diffeomorphisms in 4 space-time dimensions remains an open problem. We should 
also stress that (2.1) is not just a Brans-Dicke action in the Jordan frame 13(, subject to 
the identification of ^ with a scalar field cf>, because (2.1), as we said above, results from the 



3 



definition [3] 

s = J_ / r( (4) ^-2A c ) 

where 



G 



1*9 



pa _ ;p 



G 3 



gd 4 x-|-- 



^GL 



f = VhiKN* -h ij Nj). 



-gd 4 x, 
(2.2) 

(2.3) 



In the formula (2.2), the last two terms in the integrand are not total derivatives since G 
is variable, and hence the Euler-Lagrange equations resulting from (2.2) differ from the 
Brans-Dicke field equations. 



On relying upon the work in Ref. 10J we know that, in a spherically symmetric space- 
time, foliated by leaves E labelled by a real time parameter t, only the radial component iV r 
of the shift vector survives, and both the lapse N and iV r depend only on the (t, r) variables. 
Moreover, the three-metric of the leaves reads as [h]] 



ds 2 = A 2 (r,t)dr 2 + R 2 (r,t)(d9 2 + sin 2 



(2.4) 



where R is the curvature radius of the two-sphere r = constant, so that the space-time 
four-metric reads 

da 2 = -^N 2 - (iV r ) 2 ) dt 2 + 2N r dr dt + ds 2 . (2.5) 
Integration over 9 and in Eq. (2.1) yields therefore 



S S [R,A,G; N, N r 



L dt, 



(2.6) 



with Lagrangian (hereafter, following 10], dots and primes denote partial derivatives with 
respect to t and r, respectively) 



dr 



N~ 1 R(A - {AN r )'){R - R'N r ) - -JV _1 A(^ - R'N r ) 2 



- N[ - A- l RR" + A~ 2 RR'A' - -A~ l R' 2 + -A 



1 

2 J 



NAR 2 uNAR 2 G 2 ■ 
~2G~ Ac+ A~G^[jP~ 2 m GG + 



N r 



(2.7) 



It should be stressed that we have made a non-trivial step, i.e. the insertion of spherical- 
symmetry ansatz into the action before performing the variations that lead to the field 
equations. In general relativity, spherical reduction in the variational principle leads indeed 



4 



to the Schwarzschild solution (see comments below, in between Eqs. (2.38) and (2.42)). 
Mathematicians have realized, by now, what are the symmetry groups for which the re- 
duction fails (see, for example, the work in Ref. [jjj]). For theories with variable G, the 
rigorous proof that spherical reduction in the variational principle is admissible is an open 
problem, but we will see at the end of section 2 that the resulting Hamiltonian constraint 
is compatible with a space-time metric which, in the limit of constant G, reduces to the 
Schwarzschild metric. 

At this stage, we might write directly the Euler-Lagrange equations resulting from the 
Lagrangian (2.7). However, since the latter involves the explicit r- integration, the passage 
to Hamiltonian variables leads to a more manageable system of coupled first-order equations 
(see (2.19), (2.20) below). In the final form of the Hamilton equations one can eliminate the 
momenta and hence recover the desired Euler-Lagrange equations, if necessary, or rather 
go the other way round. One may further check the equivalence of the Hamiltonian and 
Lagrangian formalism in the relevant case of the Euler-Lagrange equation for G itself, which 
is possibly the major novelty resulting from the action (2.1). On writing the Lagrangian 
(2.7) as 

POD 

L= L(r,t)dr, 



the Euler-Lagrange equation for G, i.e. 

d dL dL_ 
Jtdd~dG~ 0, 

leads to (on setting N r = for simplicity) 

d_ ( fiAR 2 - \ _ f" 3 iikR 2 G 2 _d_ (KR 2 K C fiR^G^ 
Jt \2NG 3 )~ ~4 G 4 N 2 ~ dG \1T~G + IT^ 

But this coincides with the third Hamilton equation (2.28) (see below), by exploiting 

which is the third Hamilton equation (2.27) when iV r = 0. 

With this understanding we further remark that, by differentiating the ADM action (2.3) 
with respect to the velocities A, R and G, we obtain the momenta (cf. Ref. jlO| ) 

Pa = —N~ 1 R(R - R'N r ), (2.8) 



5 



PR 



A(R - R'N r ) + R(A - (AN r )') , (2.9) 



Po = ^^(G-G'W). (2.10) 
Equations (2.8)-(2.10) can be inverted for the velocities, i.e. 

A = -NR-\Rp R - Ap A ) + (AN r )', (2.11) 

R = -NR- 1 p A + R'N r , (2.12) 
2 NC 3 

G = ~ T1 y I PG + G'N\ (2.13) 
\i AR Z 

The ADM action (2.6) can be cast into the canonical form by the Legendre transform 

S s [A,R,G,p A ,p R ,p G ;N,N r ] = j dt j dr(p A A + p R R + p G G - NH - N r F r ). (2.14) 

The insertion of Eqs. (2.11)— (2.13) into (2.7) and (2.14) leads to two equivalent expressions 
of the Lagrangian, so that the functions multiplying lapse and shift therein must be equal. 
Hence we find (cf. Ref. (loj) 

„, v PrPa lAp A G 3 2 RR" RR'M R' 2 

1, \R\ /iG' 2 R 2 ,„„„, 
-2 A + lC A ' + i«?' (2 - 15) 

H,{r,t) = R'pg - Ap' A + G'p a . (2.16) 
In the course of deriving Eq. (2.16), we have found in p A A + p R R from Eq. (2.14) a term 

(AT)'APA = -^(N r Ap A ) - N r A'p A - N r Ap' A . (2.17) 

Thus, Eq. (2.16) holds because, upon r-integration, the first term on the right-hand side of 
q. (2.17) gives vanishing contribution subject to the fall-off conditions in Sec. IIIC of Ref. 



101 ] . The second term on the r.h.s. of Eq. (2.17) is cancelled exactly in the integrand of Eq. 



(2.14), while the third term on the r.h.s. of Eq. (2.17) leads to —Ap' A in Eq. (2.16). The 
effective Hamiltonian now reads 

/oo 
[NH + N r H r + u n tt n + z/tt;] (p, t)dp, (2. 18) 

-oo 

where ir^ and 7Tj are the primary constraints which occur because the ADM Lagrangian 
(2.7) is independent of time derivatives of lapse and shift. The general Hamilton equations 
are therefore 

j t Q(r, t) = [Q(r, t), H(t) } , Q = A, R, G, (2.19) 



d 

Jt 1 



PqM) = \Pcj(r,t),H(t) \ , p Q =pA,p R ,p G 



(2.20) 



to be solved for given initial conditions A(Q),R(Q),G(0),pa(0),Pr(0),Pg(Q) satisfying the 
constraint equations 

(2.21) 



H « 0, H r &0, 



where the weak-equality symbol ~ denotes equations which only hold on the constraint 



manifold 
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We now exploit the relations 
d 



dt 



Q(r,t) 



j t P Q (r,t) 



N(p, t) {Q(r, t), H(p, t)} + N r (p, t) {Q(r, t), H r (p, t)}j dp, (2.22) 

N(p, t) {p Q (r, t), H(p, t)} + N r (p, t) {p Q (r, t), H r (p, t)}] dp, (2.23) 
dH d , 



{Q(r,t),H d (p,t)} 



{p Q (r,t),H d (p,t)} 
and define the vector field 



9PQ 
dH d 



dQ 




(p,t)5(r,p), H d = H,H r , 
(p,t)5(r,p), H d = H,H r 



d 



dt dr 

to find, for all values of lapse and shift, the general Hamilton equations 

D tr A « NU U D tr R « NU 2 , D tr G « NU 3 , 



having set 



Vi = 



Vo = 



Pa 



G 3 



2R 2 pcA 2 R 
PrPa , Api . 2G 3 



D trPA « NV U D trPR « W 2 , ArPG « NV 3 



p R ApA 
rr _ PA 

U2 = -R> 
_ 2G 3 

^ = M^ PG ' 

d_fRR^_RRfA^ R?\ 1 _ //G' 2 fl 2 
3A V X A^ + 2A/ + 2~2G C+ 4A 2 G 3: 



2^G 



+ 



i? 2 i? 3 ^Ai? 3 



9 fRR" _RR'A' R' 2 \_AR _ /iG' 2 R 
dR\~A A^ + 2Ay~G" c ~2 AG 3 



V* 



3G 2 



fiAR 2PG ~ dG [ ~1T~G + 4 T G 3 " 



(2.24) 
(2.25) 

(2.26) 

(2.27) 
(2.28) 

(2.29) 
(2.30) 

(2.31) 
(2.32) 
(2.33) 
(2.34) 
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The six equations (2.27) and (2.28) should be studied, for given initial conditions, jointly 
with the two constraint equations (2.21) and with the ADM relations for lapse and shift in 



the (£, r, 9, 0) coordinates, i.e. [lGj 



A = V-FT' 2 + F^R' 2 , (2.35) 
A7 R'f-T'R 

N = , (2.36) 

AT = —FT'f + F~ l R'R 
A 2 

It is reassuring to note that such equations make it possible to recover the Schwarzschild 
solution in general relativity. For this purpose, it is enough to choose the foliation defined 
by 

R(r,t) = r, T(r,t) = t, (2.38) 

for which A^ r = and hence p^ pa 0, pn ~ 0, with Hamiltonian constraint H pa reducing 
to (from Eq. (2.15)) 

RR A' 1 . i o \ / 

— ~2A (1 - A) - (2 - 39) 
Such an equation is solved by R = r and 

A = F- 1/2 = M — J . (2.40) 

Moreover, the weak equation (2.39) can then be used to cast Eq. (2.32) in the form 

2RR'M RR' OA' R' 2 1 

V\ pa 1 1 h - 

1 A A 2 A 2 dA 2A 2 2 

K -!(2X-i) + 2F^ A - A, ) + ^ + ^°- p- 41 » 

Along the same lines, we obtain the weak equation 



V, pa -— 



d (RR" RR'M R' 2 



OR V A A 2 2A 
AdR dR K "" ' " 2^ " ' ' 2 



(RR')-R' 2 -l(l-A 2 ) + h> ' 1 



--X^^ + 2l^>°- ^ 
Our remark agrees with the findings in [if]], but our analysis offers the advantage of not 
having to eliminate R and from the Hamiltonian analysis, which is important when G 
and A c are allowed to vary. 
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In the latter case, we can solve the general Hamilton equations expressed by (2.27)-(2.34) 
when the space-time foliation is again given by (2.38) with vanishing shift and p\ ~ 0, pn ~ 
0; Pg ~ 0. In this static case, where only the spatial gradient of G is nonvanishing, we 
assume a fixed-point relation as in [zj]: 

A C G = A c (r)G(r) = constant = k, (2.43) 

so that the Hamiltonian constraint H pa yields now 



RR'M 



—-(1 - A 2 ) 

2A V ; 



kAR 2 fi G' 2 R 



(2 44) 

A 2 ' 2A y ~ " ' ' 2G 2 ' 4 AG 3 ' K ' ' 

The functions V%, V2 and V3 defined in (2.32)-(2.34) are again weakly vanishing, since, by 
virtue of (2.38), (2.43) and (2.44), 



Vx 



+ 



2 

~A 
]_d_ 

AJdA 
R' 2 



1 

2A 1 



kAR 2 ^ G' 2 R 2 
2G 2 + 4 AG 3 



A 2 ) + 



'^(1-A 2 
2 v 



fcA 3 ^ 2 uAG' 2 R 2 



2G 2 



4 G 3 



1 

2A 2 + 2 



fc/g 2 /i G /2 fi 2 
2G 2 " + 4 A 2 G 3 



(2.45) 



OR 
0, 



fcAi? 2 /i G" 2 i? 2 



2A 



(1-A 2 ) + 



+ 



2G 2 4 AG 3 



kAR fi G' 2 R 



G 2 



2 AG 3 



[2.46) 



V* 



d_ 
dG 
_d_ 

dG 



AR 2 A c fi R 2 G' 2 
~^T~G + IT G 3 " 
RRA' 



A 2 



A 2 



0. 



(2.47) 



If no infrared fixed-point relation such as (2.43) can be assumed, Eqs. (2.44)-(2.47) still 
hold provided that one replaces k therein by the product A C G. For example, the Hamiltonian 
constraint (2.44) takes the form 



RR'M 



AR 



M+^7 



fiG ,2 R 2 



(2.48) 



A 2 2A y ' 2G 4 AG 3 
Moreover, it is always true that A c and G are not independent variables but are function- 
ally related. This is clearly proved in the Hamiltonian framework advocated in our paper. 
Suppose in fact that A c were an independent dynamical variable. The primary constraint 
of vanishing momentum conjugate to A c would then be preserved in time, from (2.1), pro- 
vided that either the lapse function or the determinant of the induced three-metric vanishes, 
leading therefore to a complete 'collapse' of the ADM geometry. 



III. NONLINEAR DIFFERENTIAL EQUATION FOR THE NEWTON PARAM- 
ETER 



We now also assume that the departure from general relativity is not so severe, so that the 
A function keeps its functional dependence on G(r), at least approximately. More precisely, 
since A should reduce to (2.40) in the case of constant G, it can only differ from (2.40) by 
terms involving explicitly the gradient of G. A good 'a posteriori' check of any approximate 
solution for G is therefore whether it has a gradient with negligible effects. We thus insert 
into Eq. (2.44) the form (2.40) of the A function, with G taken to depend on r only, and we 
obtain eventually the nonlinear differential equation 

A{r)G' 2 {r) + G 2 {r)G'(r) - B(r) = 0, (3.1) 

where (hereafter, G n is the Newton parameter on solar system scale) 

A(r) = ^r(l- ^-rl G n < 0, (3.2) 



2 V. 2MG(r) 
B( r ) = —r 2 c 2 G n > 0, (3.3) 

and we have restored the physical units, since we want to make estimates on real situations. 
The statement A < is true only for /i > and r sufficienly large, which is surely true for 
normal astrophysical objects, like sun and galaxies. The case of the immediate neighbour- 
hood of a black hole horizon is more involved and goes beyond the aims of the present paper 



(cf m- 

We therefore see that Eq. ( 13.1 ft admits always two distinct positive values for G'(r). This 
means that in any case G(r) is monotonically increasing with r. We shall see in a moment 
that this can be reconciled with the requirement that the metric should be of Minkowski 
type at infinity. Another problem is posed by the two disjoint solutions. We assume however 
that fi is sufficiently small to get the first term of Eq. (13.11) negligible in our case. We have 
thus to treat a much simpler equation (a more accurate treatment and justification of this 
assumption is given in the next section), i.e. 

G 2 dG = ^-%^r 2 dr, (3.4) 

which leads to the growth of the Newton parameter according to 

GW=G„(l + ^)', (3.5) 
10 



where we have set the integration constant, in order to obtain a correction to Newton's law, 
as desired. 

At large r, the Newton parameter obeys therefore the approximately linear relation 

Therefore, reverting to Eq. (I2.37p . we see that the function F = (1 — 2GM/rc 2 ) tends 
asymptotically to the constant 

F=1 .(^y, 0.7) 

so that, by rescaling appropriately distance and time, we may obtain flat space. 

A very interesting feature of Eq. (I3.6P is that it gives just the correction necessary to 
obtain perfectly flat rotation curves of galaxies. Let us indeed rewrite it as 

G = G n (l + ar 3 g M; 1 p 3 g ) 1/ \ (3.8) 

where M g is the galaxy mass, a = and p g = r/r g is the distance rescaled according with 
a typical galaxy length r g . We see that the correction is effective at say p g ~ 1.2, if we take 
aTgMg 1 ~ 0.5. At the solar system scale we get (p s = r/r s ~ 1) 

G = G n (1 + arlM; l p^) 1 " = G n (l + ' ^ G n (l + 10- w p 3 s f 3 , (3.9) 

and the correction is absolutely irrelevant. 

It is also interesting to note that, at large r, we get for the radial velocities of galaxies 
the relation 

^^MocM 1 / 3 , (3.10) 

where the proportionality constant is equal to (fcc 2 G ? n /2) 1 / 6 . Now, the usual theoretical ex- 
pression for the Tully-Fisher relation is v oc M 1 / 4 , and is computed with the usual Keplerian 
law for velocities. On the other hand, we have the observational relation 

Mag = -7.681og 10 ( -^-v ) - log 10 g, (3.11) 
\sini / 

where Mag is the absolute magnitude, i is the visual angle of the galaxy and q is a number 
which depends on the optical band 18]. The first coefficient, which is the only relevant 



11 



one for our purpose, has small dependence on the band. If we consider also the theoretical 
definition of magnitude, i.e. 

Mag = -2.51og 10 L + log 10 a, (3.12) 

where L is the absolute luminosity, assumed proportional to the mass, and a is again de- 
pendent on the band, but irrelevant, we obtain eventually 

v oc M 2 - 5 / 7 - 68 = M 0325 ~ M 1 / 3 , (3.13) 

n 

so that we obtain a striking agreement with our Eq. (3.10), unlike the work in [5], where 
no agreement with the empirical Tully-Fisher relation is found (see comments in the last 
paragraph of section 4 therein). The reason for this improvement lies in the fact that, in 
our case, the correction (3.5) to the Newtonian G is not parametrized only by universal 
constants, but by the mass of the gravitating source of the field. What is instead depending 
only on universal constants is the proportionality coefficient v/M 1 ^ 3 in (3.10). 



IV. QUALITATIVE ANALYSIS OF THE EQUATION FOR G 

Let us now consider again Eq. (13. ip and show that indeed its replacement with the much 
simpler Eq. (13.41) is justified. First, let us point out that, for any astrophysical object 
different from a black hole we may write safely 

* r > s -mw, G - (4 ' 1) 

Then Eq. ( 13. II) may be rewritten, in c = 1 units, as 



= G 3 (r)± > /G«(r)-2a^G(r) 
ar 2 

where a = fiG^/rg, (3 = kG 2 n jr s ^ and r s = 2G n M is the Schwarzschild radius of the object. 

We are considering only one of the two equations generated by Eq. (13.11) . the study of the 
other being made along the same lines. We see that we may reduce the number of relevant 
parameters to only two. This equation can be exactly solved in the cases a = and (3 = 0. 
The first corresponds to the solution examined in the previous section, which we rewrite 
here as 

G(r) =/3 1/3 (r 3 + 3w) 1/3 , (4.3) 
12 



where w is an integration constant. 

The second case is possibly even more interesting, since it corresponds to chosing A c = 0, 
which is closer to the Schwarzschild geometry. There is then a trivial solution G = constant, 
as well as 

G(r) = ^ IT ^- at r»r s , (4.4) 

where x is the integration constant. We see that this solution also tends asymptotically to a 
constant, and hence satisfies the consistency check stated at the beginnin of section 3. If this 
regime is reached sufficiently late, an emulation of dark matter might be obtained again on 
taking a linearized approximation of (4.4) (cf. comments in section 5). Everything depends 
on the values of the parameters and confrontation with observations. 

Let us now show that the first solution dominates at large r. For this purpose, let us 
consider the ratio of the first to the last term in Eq. (13. II) . i.e. 

6 = wWy (4 ' 5) 

It is clear that our approximate treatment will be good as long as 5 << 1. We substitute in 
S the approximate solution and obtain 

s >= a ^' 3 w^w - (4 ' 6) 

The behaviour of this function is independent of a and (3. It tends asymptotically to zero 
and has a maximum at r m ~ 2.29i/; 1//3 . Therefore, provided we start the integration at 
r start > r m, if the approximation is valid there, it will be increasingly accurate as r gets 
larger. 

On the other hand, let us suppose that at r start the opposite occurs, and <5(r start ) >> 1. 
We may thus substitute the other solution, finding 

= 2a 3 / 2 

2 /3r 3 / 2 (4 + 2xr) 5 / 2 ' { ' ' 

We see that again the behaviour of 62 is independent of a and (3, but (which is most 
important) we always have that 8' 2 < 0. Therefore, even if at the beginning S(r start ) >> 1, 
as r increases the condition is reversed and we may say (approximately) that, when <5 2 = 1, 
we may switch off this solution and revert to the first one, which prevails asymptotically. 
The intermediate situation is of course somewhat delicate, but a numerical analysis, made on 
the full equation, with suitable choice of the parameters involved, confirms these statements. 
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V. CONCLUDING REMARKS 



In the first part of our paper, we have extended the Hamiltonian analysis of spherically 
symmetric gravity [h], to the case of variable Newton parameter and variable cosmolog- 
ical term, obtaining eventually the non-linear differential equation (3.1) for G(r), under the 
non-trivial assumption that Eq. (2.40) can be taken to hold. We have then shown that the 
treatment of A c and G as dynamical variables, together with the fixed-point condition, gives 
encouraging chances of emulating the presence of dark matter in long-range gravitational 
interactions, at least at galactic scale. Several open problems should be now studied, i.e. 

(i) The legitimacy of the fixed-point assumption. 

(ii) The validity of Eq. (2.40) when G depends on r. 

(iii) The detailed numerical proof that also our solution (4.4) with vanishing cosmological 
constant can fit the flat rotation curves of galaxies. 

(iv) Can our Hamiltonian approach make it possible to study weak-lensing observations, 
ihat are recently found to provide a direct empirical proof of the existence of dark matter? 



(v) How to perform the Hamiltonian analysis with variable G in the small-r region which is 
relevant for black-hole physics? 

The main source of future developments is the confrontation of our theoretical results with 
observational data. This is not a simple task, because a galaxy is not a pointlike source. 
Therefore a separate paper is in order on this topic as well as the other open problems listed 
above. Anyway, since the free parameters of our theory, a and /3, depend on \x and k, which 
are completely undetermined at the moment, we hope to be able to obtain reliable numerical 
solutions (see Figs. 1 and 2), appropriate for comparison with observations. 
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FIG. 1: The two approximate solutions of Eq. (4.2) and a rough attempt to join them appropriately. 
The red curve is obtained in the case a = 0, the green is for = 0, so that the red is appropriate 
when 5 « 1 and the green when 5 » 1. The black dashed is obtained by using the first solution 
for 5 < 1, which happens at the beginning and at the end, and the second in the other case. The 
discontinuity in the derivative results from the rough procedure. This curve should be compared 
with the ones of Fig. 2, obtained numerically from the full equation. 
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FIG. 2: Numerical integration of the full Eq. (4.2) when a = 0.1,(3 = l,r sta rt = 1; for various 
values of G s tart- Units are arbitrary. Comparison with Fig. 1 shows that the qualitative analysis 
is correct. 
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